arXiv:hep-ph/9808447vl 27 Aug 1998 


IASSN S-HEP-98/74 
OSU-HEP-98-5 
lhcp-ph/ 9808447 


A Reexamination of Proton Decay 
in Supersymmetric Grand Unified Theories 


K.S. Babu 1,2 and Matthew J. Strassler 1 

1 School of Natural Sciences, Institute for Advanced Study, Olden Lane, Princeton, NJ 08540, 

USA 

2 Department of Physics, Oklahoma State University, Stillwater, OK 74078, USA.* 

babu@osuunx.ucc.okstate.edu, strasslr@ias.edu 

Abstract 


We reconsider dimension-five proton decay operators, making semi- 
quantitative remarks which apply to a large class of supersymmetric GUTs 
in which the short-distance operators are correlated with the fermion Yukawa 
couplings. In these GUTs, which include minimal 517(5), the operators 
(it c d£)t (djU T ), induced by charged Higgsino dressing, completely dominate 
for moderate to large tan/3. The rate for p —> (K + ,tt + )17 t grows rapidly, as 
tan 4 /?, and the K + to ir + branching ratio can often be precisely predicted. 
At small tan (3 the operators (udi){djv) are dominant, while the operators 
( udi){ul~ ), with left-handed charged leptons, are comparable to the neutrino 
operators in the generic GUT and suppressed in minimal GUTs. Charged- 
lepton branching fractions are always small at large tan f3. The electron to 
muon ratio is small in minimal GUTs but can be larger, even of order one, 
in other models. All other operators are very small. At small tan /? in non- 
minimal GUTs, gluino and neutralino dressing effects on neutrino rates are 
not negligible. 
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Unification of the three gauge forces of the standard model into a single symmetry group 
is an attractive idea with a long history UIJ. Grand Unified Theories (GUTs) provide a 
beautiful explanation of the multiplet structure of the standard model fermions, and predict 
the value of one of the gauge couplings at the weak scale in terms of the other two. The 
success of this prediction in the context of the minimal supersymmetric extension of the 
standard model (MSSM) has given new impetus to the study of SUSY GUTs over the past 
decade; for a review see 0 . 

GUTs place quarks, leptons and their antiparticles in the same multiplet, thereby mak¬ 
ing the nucleon unstable. The momentum scale of grand unification, as inferred from 
the extrapolation of the standard model gauge couplings, is M GUT ~ 2 x 10 16 GeV, so 
large that nucleon decay rates are compatible with present experimental limits. In su¬ 
persymmetric GUTs there are two sources of baryon number violation. Dimension-six 
(four-fermion) operators are generated by the exchange of superheavy gauge bosons of the 
GUT symmetry group. With M G ut ~ 2 x 10 16 GeV, the expected proton lifetime from 
these operators is r(p —> 7i°e + ) ~ 10 36±2 yr. 0# the present experimental lower limit is 
r(p 7r°e + ) > 2 x 10 33 yr. [|1(J. The second source of proton decay — the subject of this pa¬ 
per — is the dimension-five (two-fermion/two-sfermion) operators induced by the exchange 
of color-triplet Higgsinos, the GUT partners of the MSSM Higgs(ino) doublets | pT| . pU|| . The 
associated amplitudes scale as Mq ut , but are suppressed by light fermion Yukawa couplings. 
They are also suppressed by a loop factor, which arises from the “dressing” of the operator 
by a gaugino or Higgsino that converts the two sfermions into light fermions. The size of the 
loop integral is highly uncertain, due to the unknown masses of the supersymmetric parti¬ 
cles. Additional uncertainty stems from the unknown values of tan (5 = (H u ) / (Hj) and the 
triplet Higgsino mass. (In branching ratios, some of the unknown parameters cancel, so they 
can be more reliably predicted than the overall rate.) Since the strength of these operators is 
proportional to Yukawa couplings, the dominant modes are those with strange quarks - the 
heaviest fermion into which the proton can decay kinematically. In minimal SUSY SU( 5), 
including all the uncertainties, the lifetime may be estimated as r{p —> K + u) = (10 28 —10 34 ) 
yr. This range is nearly eliminated by the experimental limit r(p —> K + i 7) > 5 x 10 32 yr. Ju|], 
suggesting that proton decay must be discovered imminently if the idea of grand unification 
is on the right track. 

These issues have been studied extensively in minimal SU( 5) 

00] for more references to the early literature — as well as in other GUTs 
calculations have been specific to particular models. In this letter we attempt to make 
statements which, while more qualitative than quantitative, apply to a wide class of theories. 
Our purpose in this paper is to summarize and clarify the expectations concerning branching 
ratios in a large ensemble of supersymmetric GUTs, in which the dimension-five operators 
are roughly correlated with the fermion Yukawa couplings. We will make this more definite 
below. Our approach is somewhat similar in spirit to that of |[V] . 

General statements of this type require a thorough analysis. To this end a systematic 
method has been developed in which rough upper bounds on baryon-number-violating op¬ 
erators may be established; this will be presented elsewhere |27j . These bounds apply under 
the following conditions. The dimension-five operators take the form 


F4H20H 


1M1- 


- see 
Most 


— 


d 2 9 {A ijkr QiQ 3 Q k L r + A ijkr U*D c r U*E c k } + h.c. 


( 1 ) 


2 

















(Capital letters denote superfields; Ui contains the squark u t and the quark u j, etc.) We 
take as our ansatz the relations Xijkr ~ \jkr ~ jf e i e j e ki independent of the index r, where 
M is an overall mass scale independent of i,j,k,r and e 3 ~ 1, e 2 = e ~ 1/25 and ~ e|. 
(The motivation |28|j29|] for the ansatz is associated with the fact that, near the GUT scale, 
most of the flavor structure of the theory can be roughly encoded in powers of e |I0|.) In 
i) the Xijkr are expressed in the gauge basis, but it will be more convenient to work in 
the “supersymmetric basis”, in which the matter fermions are expressed as mass eigenstates 
and the sfermions are expressed as the supersymmetric partners of those eigenstates. To 
convert from the gauge (g) to the supersymmetric (m) basis, we must rotate the superfields 
as U- fl} = R^Uj m \ D\ g ' = RfjDj m \ etc. so that the fermion mass matrices are diagonal. In 
the supersymmetric basis, 


W dim -s = C>t jk r U ?DjUlE r - KjkrUi D'jDINr + % iikr Uf D?U]E k )e a/h . (2) 

where a,/3 ,7 are color indices, we have dropped the superscript (m), and we have defined 
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the coefficients X'- ]kr and Xijkr are defined analogously. As part of our ansatz, we assume 
the matrices R u , R d , R u , R e — those associated with the 10 representations of SU( 5) - 
have the same texture as Vckm ; we do not assume this for the other R matrices. If Vckm 
satisfied ~ min {d/ej, Cj/ei}, then we would have X\ jkr ~ X^ kr ~ X ljkr ~ X ljkr ~ e^jCk- 
However, the Cabibbo angle 6 C is of order 5 ei/e 2 and can enter into the A. An enhancement 
by a factor of 0 c /(f-) ~ 5 may occur for each index i,j, k taking value 1 |50); thus 
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Note Sn + + Ski < 2 because of the antisymmetry properties of the A coefficients. More 

details of this ansatz are given in J77J and 0 - 

If the coefficients A in a given GUT are all of order or less than those appearing in Eq. (|4|), 
then we can apply the results of (27J to the model. Although we will not classify them here, 
many GUTs (including minimal 517(5) and 50(10) and their more realistic variants, as well 
as the ten-centered models of |28 H30|1 ) are compatible with the ansatz (||). We will now 
demonstrate that the A coefficients in the minimal 517(5) GUTs are consistent with (|f). 
In the supersymmetric basis, the superpotential induced by the color-triplet Higgsinos in 
minimal 517(5) is 

ll'r = (i UiDje^fiVij + UfEff^Hc + (DXR, - U.E^hj + U^e V~h ■)IIc . (5) 

Here He and He are the color triplet and anti-triplet from the 5 H and 5 h of the Higgs 
multiplet, V is the CKM matrix, /, (hi) is the diagonal Yukawa coupling matrix of the up- 
quarks (down-quarks and leptons), and the cq are phase factors with X^ =1 T: = 0. Exchange 
of the Higgs color triplets, of mass M, leads to the dimension-five operators 
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along with other terms irrelevant for proton decay. 

If it were the case that 9 C ~ ei/e 2 , so that V tJ ~ min{ej/ej, e^/e*}, along with f t ~ ef, 
hj ~ €i(, where ( = tan/3/60, we would have 


fiVijSkrhr ~ e*€r5 kr mm fa/ 6j , €j / 6i}( < 6^6^ 
fiVijV kr h r ~ ele r mm{ei/eji6jl6i}mm{e k /e r ,e r /e k }C, ^ e^-CfeC 
fi V jkV*h r ~ Cjfr min{ej/e fe , e fc /e,,} min{e;/e r , e r /e0C ^ (7) 


All coefficients in (|jj) would then be equal to or smaller than those in Eq. (£|); note the 
overall factor of ( can be absorbed into the overall mass scale which affects all A equally. 
Now let us account for the fact that 9 C ~ 5ei/e 2 . Inspection of Eq. (|7|) shows that factors of 
the Cabibbo angle can only enhance a coupling above e^e*, if one of the indices i,j, k takes 
value 1 ; r = 1 cannot cause such an enhancement because of the factor hi, which reduces 
, Ajjfci below its ansatz value. If two such indices take value 1, then one may get a 
double enhancement. The resulting coefficients are the same as or less than those appearing 
in Eq. 0), and thus the upper bounds of |27j] are applicable. 

Rough upper bounds from [^7] on various four-fermion baryon-number-violating opera¬ 
tors are given in the Table. All operators not shown are negligibly small. The Table shows 
the different contributions from various gauginos and Higgsinos. Wino-Higgsino mixing is 
not listed; in each case the upper bound on such contributions lies between or below the 
pure Wino and the pure Higgsino bounds. Explicitly indicated are factors of e, 9 C , gauge 
coupling constants g* and third-generation Yukawa couplings y bl y T which stem from left- 
right mixing or Higgsino couplings; note y b and y T are roughly of order (. Also appearing 
are parameters 7 and 5. The first measures the extent to which the b, t are split in mass 
from and mixed with the other left-handed squarks; if the messenger scale of supersymmetry 
breaking is high, then 7 ~ 1, while if it is near the weak scale 7 may be small. The factor 
5 = (Acos/3 + fi sin/3)(A sin/3 + //cos (3)v/fh 2 parameterizes the size of left-right squark 
mixing; here v, p, m, A are 246 GeV, the Higgsino mass, the universal squark mass and the 
coefficient of the trilinear scalar terms. An overall factor of M' 1 // is omitted from every 
entry. Neutrino flavors have not been distinguished, and the flavor of a charged lepton is 
only indicated when the bounds depend on the lepton flavor. For additional details see [|27| . 
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Rule-of-thumb upper bounds on coefficients of four-fermion baryon-number-violating op¬ 
erators, as found in using the ansatz in Eq. (ff\). Important loop factors, matrix elements 
and overall coefficients are omitted and must be accounted for when using this Table. See 
the text for further explanation. 

At small tan/3 the operators udsv , uusp and uuse have the largest bounds, while at 
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large tan/3 the operator ( u c d c )^(sv) 


potentially exceeds all others. The replacement of the 

As we will see, it is always true 


s quark by a d quark engenders a small suppression l|14]Jl5| 
that many of these operators saturate or nearly saturate their bounds. It follows from the 
Table that operators with right-handed leptons, even u^p^su, are negligible at all values of 
tan (3, as are operators with a left-handed charged lepton and two right-handed quarks. We 
will therefore restrict out attention to the first four rows in the Table. 

At small tan f3 the largest four-fermion operator is udsv , induced through dressing of the 
CDSAfTD /iAy, term in the superpotential by the charged Wino, as in Fig. |l] [[To| . p~8| . In 
minimal GUTs the second-generation amplitude is proportional to y c y s d 2 ( y c is the charm 
quark Yukawa coupling, etc.), while the third-generation effect is roughly of the same order. 
Since y s y c 9 2 c ~ ~ e 3 # 2 ^, the coefficient of udsv saturates its upper bound in 

the Table. Similar statements apply to uddu , with an additional factor of 9 C . Interference 
effects between second- and third-generation diagrams [which unfortunately depend on the 
otherwise-unmeasurable phases <7,; in Eq. ([])] might suppress either p —> K + D or p —> tt + u 
l8| but not both ||. Accounting for this and for the different hadronic matrix elements 
and phase space factors, one folds the branching ratio T(p —> K + v)/T(p —> n + T7) tends to 
lie between .1 and 1. (Throughout this letter we assume that sfermions of the same charge 
are degenerate at the messenger scale of supersymmetry breaking.) In more general GUTs 
satisfying Eq. (||), even this predictivity may be lost, though the tendency is still the same. 
The rates for these processes grow as tan 2 (3. 

At large tan /3, the largest operator in these GUTs is [u c d c y (sis) , first discussed in |p~0l . [24 


This operator is obtained through Higgsino exchange or through left-right sfermion mixing, 
which always entails a factor of the associated right-handed fermion mass. The largest 
effect therefore comes from changing a third-generation right-handed squark or slepton to 
a light left-handed quark or lepton. Specifically, begin with the operator U c D c T c t c in the 
superpotential, and then convert the sfermions t c f c to the fermions sz/ T . The Higgsino and 
Wino dressing diagrams are given in Fig. |2|; there are also diagrams with Higgsino-Wino 
mixing whose size is intermediate between them. The diagrams are proportional to the 
coefficient A 1331 . For the Higgsino dressing one has the couplings ytV t * ~ y t e 2 for the t c — s 
vertex, y T for the f c — v T vertex. A similar structure emerges for the Wino diagram but 
suppressed by gauge couplings and mixing factors; the Higgsino diagram usually dominates. 
The coefficient of this operator in minimal SU( 5) is proportional to //d// 2 // T VL; ~ e 3 ([ 2 . This 
lies e/9 c below its upper bound, because, as is also the case in many other GUTs, A 1331 = ydyt 
is a factor of e/9 c below the ansatz in Eq. ([|). This will not be true in non-minimal GUTs 
where the matrix R d is non-hierarchical and i?“ 2 ~ 9 r , as in certain ten-centered models flSH 


where the bound on this coefficient can be saturated. 

The same set of graphs, with only the external down-type quark flavors changed, gives 
{u c d°y (su T ), ( u c s c )\dv T ) and ( u c d c )\dv T ). In minimal SU( 5 ) the three amplitudes are in 
the ratio ydV ts , ysV us Vtd, ydVtd', note these ratios are independent of the phases cq, unlike the 
udidjU case. The coefficients of {u c d c y (su T ) and (u c s c y (dv T ) are of opposite sign, and may 
be of the same order. (This is slightly inconsistent with the Table; however, the estimates 
therein are rough, and, as noted above, the (u c d c y (sv T ) amplitude does not saturate its 
bound.) However, the hadronic matrix elements of these operators have opposite signs (to 
see this use ||Ti] . |T^| ) so interference in p —» K + T7 r is constructive; also, since {u c d c y (sv T ) has 
a significantly larger matrix element, its contribution dominates. The operator {u c d c y (du T ) 
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leads to p —> i t + v t ] since its coefficient is in a known relationship to the other two, a precise 
prediction for the branching ratio F(p —> 7 r + I 7 T )/r(p —> K + i 7 r ) is possible, independent of 
the supersymmetric spectrum. This predictivity is retained in models where A 1331 /A 1332 is 
determined. Even this information is unnecessary in GUTs where A 1331 is as large as allowed 
in Eq. (|j), since in this case the (u c s c Y (dv T ) operator is negligible and only A 1331 appears in 
the amplitudes. 

Since these amplitudes go as ( 2 , the rates for p — > K + v T , tt + v t grow like tan 4 /?, leading 
to short lifetimes and corresponding strong constraints at large tan /?. We fold that the 
amplitudes for udsv and ( u c d c )\sv ) become comparable in most GUTs for tan/? somewhere 
between 3 and 15; in a minimal supergravity SU( 5) GUT (accounting for the different short- 
distance renormalizations and hadronic matrix elements of the two operators) we fold this 
number is of order 9(m^/p,), with large uncertainties from the ratio yd/ Vs, the sfermion 
spectrum, poorly measured CKM angles, and the third-generation contribution to the udsv 
operator. [] For tan /? ~ 60, the rate from ( u c d c y(sv T ) dominates that from udsv by 10 — 1000. 
In those non-minimal GUTs where the bound on (u c d c Y (sv T ) is saturated, the amplitude 
is enhanced by another factor of 9 C and can dominate for even smaller values of tan /?. The 


constraints on large tan /? models from these effects, although discussed in |24j , do not appear 
to have been fully incorporated in the literature. 

The bounds on uusi and uudi are comparable to those on udsv and uddv. However, in 
minimal (and some non-minimal) GUTs, these bounds are not saturated. From the Table, 
we see that at small tan/3 we need only consider Wino exchange. The argument that uusi 
is highly suppressed [[H|,[nj is that all contributions are proportional to /1 = y u ~ e 4 , which 
makes the resulting amplitudes of order f\h 2 ~ e 5 , smaller than the e 3 9 2 bound. To see this, 
consider the contribution of UiDjDkAT r in Eq. (H); the sneutrino must couple to the Wino, 
implying i = 1 and giving a factor of f\. If instead we use UiDjUkE r , either i — 1 or k = 1. 
The former case gives /1 directly, and the latter gives /1 through a unitarity cancellation: 
the diagram in Fig. El is proportional to 


3 

'12 E 

i,j =1 


W 2 >: = h 2 uv* 2 v lp . 


( 8 ) 


However, the unitarity cancellation in Eq. (||) partly fails due to subtle renormalization 
group effects. The operator UiDjUj.E r is proportional to V t3 at the GUT scale, but after 
renormalization to low-energy it is no longer proportional to V %3 at the weak scale. This 
effect is of order (1 — y 2 /y 2 ) 1 / 24 p4| , |35| (here /// ~ 1.1 is the fixed-point value of y t ) or about 
0.1 for tan/? ~ 1.4 — 3. Specifically, consider the UiDjHc coupling F i3 (we use the notation 


of p 8 || ;) naively F %3 = but in fact F 31 , F 32 differ from this at the weak scale by ~ 10%. 
The sum over j in (||) becomes h 2 Vi 2 J2j Uq F tJ V ip ~ (0.1)h 2 Ui 2 % 3 % 33 %ip .f'i , which is about a 
tenth of h 2 j\ . The loop factor for this effect can be enhanced if the third-generation squarks 
arc lighter than those of the first two generations; still, even if the CKM angles, /1 and the 


*While this letter was in preparation, a preprint appeared which studies the (n c d/)f (djv) operators 
in the minimal SU(5) GUT [33|. 
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spectrum are at the edge of their ranges, it can never be as large as the leading contribution. 
A second effect of the same order (if the messenger scale of supersymmetry breaking scale 
is high) comes from the mixing and mass splittings between the down-type squarks. Both 
of these effects, while interesting, are most likely lost in the uncertainties surrounding y u . 

In many non-minimal GUTs, the unitarity cancellation in (| 8 ]) simply does not occur. A 
sufficient condition for this is that F^ be hierarchical but not precisely equal to fiVij. If Fij ~ 
/i min{ej/ej, Cj/ej} as occurs in many realistic models of flavor (including those in which 
higher-dimension operators or non-minimal Higgs bosons contribute to the up-type quark 
masses), the sum VyF tJ does not equal fidu, instead it gives d c f 2 ,ed c f 3 for i = 2,3. In this 
case the diagrams in Fig. |3] with i — 2,3 are proportional to cf'^V-'n-, h 2 V^ 2 {e 9 c f 3 )V 3 2 ~ 

e 3 6 * 2 . This saturates the bound in the Table, and thus in these theories the neutrino to 
charged-lepton branching ratio is order one at small tan/3 (although the hadronic matrix 
elements favor neutrinos.) Enhancements of F i3 by factors of do not change this 

conclusion. 

The charged lepton rates from these operators increase only as tan 2 (3. As suggested in 
the Table, at large tan /3 a bigger contribution to these operators may come from up-type 
squark mixing in gluino dressing pU| . |?Tj| ; however the branching fraction to charged leptons 
remains small due to the large (u c d°y su T operator. 

It is clear from the Table that all observed muons should be left-handed, in contrast to 
proton decays mediated by dimension-six operators |36|,|37|. Although the Table suggests 


branching fractions to electrons can be of the same order as those to muons, in minimal and 
many non-minimal GUTs they are much smaller |T4| , |T5|| . In such GUTs the coefficients A|- fcr 
are not roughly independent of r, in contrast to the ansatz (£|). However, the suppression 
factor is model dependent, and there are theories (such as ten-centered models |28 H30| 1 ) in 
which electron and muon decays do have comparable rates. The electron-to-muon branching 
ratios thus are good probes of flavor physics ] 2 B| . 

Finally, as evident in the Table, dressings involving gluinos and neutralinos can be im¬ 
portant at small tan/3 in their contribution to the udsv,uddv operators. In minimal GUTs 
these contributions are suppressed by y u , but in non-minimal GUTs (as in ||24[) they may 


become important. Gluino dressing is naively subleading due to the symmetry structure of 
the dimension-five operators (b||T(|,|T7j , and can only play a role when there is significant 
flavor violation [^0].|TJ . In the supersymmetric basis, flavor violation appears as intergener- 
ational squark mixing. At small values of tan/3, d* — dj mixings are induced proportional to 


VtVxV$ 


~ ei6j 


times the factor 7 in the Table. I 11 minimal SU(5), because neutral gaug- 


inos do not change flavor and intergenerational mixing is small in the up-squark sector if 
tan/3 <C 60, the udsv and uddu operators must come from UiDjD^Mr with 3 = 1; from 

giving effects proportional to y u y s ~ e 5 -C e 3 # 2 . 


~ 6 


13 y and both R ™2 and R 21 are of order 9 C , then it is 


Eq. (||) this implies a factor of f\ = y, 

But in a non-minimal GUT, if F t 
possible that Fn ~ / 26 12 , giving effects of order e 3 # 2 , comparable to the leading Wino dress¬ 
ing contributions (24J. Dressing by neutralinos (HVino and H-ino) is similar to the gluino 
dressing; however, squark mixing is not required for neutralinos to contribute to udsu, so 
their effects may overshadow the gluino if 7 is small [^7| . Gluino dressing also can contribute 


without mixing to udsv if there is substantial D-term splitting of squark masses 
Further details will appear in [^7|. 
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FIGURES 



FIG. 1. Wino-dressing diagrams contributing to udsv\ there are also diagrams with s,d ex¬ 
changed. 



FIG. 2. The Higgsino- and Wino-dressing diagrams leading to u c 'd c 'sv T . 



FIG. 3. The important diagrams contributing to uusfi. 
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